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Abstract 

Random Hermitian matrices with a source term arise, for instance, in the study of non-intersecting 
Brownian walkers [1, 20] and sample covariance matrices [4]. We consider the case when the n x n 
external source matrix has two distinct real eigenvalues: a with multiplicity r and zero with multiplicity 
n — r. The source is small in the sense that r is finite or r = C(n 7 ), for < 7 < 1. For a Gaussian 
potential, Peche [29] showed that for \a\ sufficiently small (the subcritical regime) the external source 
has no leading-order effect on the eigenvalues, while for \a\ sufficiently large (the supercritical regime) 
r eigenvalues exit the bulk of the spectrum and behave as the eigenvalues of r x r Gaussian unitary 
ensemble (GUE). We establish the universality of these results for a general class of analytic potentials 
in the supercritical and subcritical regimes. 

Contents 



1 Introduction 2 

1.1 Definition of the supercritical, subcritical, and critical regimes 3 

1.2 The kernel and its connection to multiple orthogonal polynomials 7 

1.3 Assumptions on A and V(z) and results 8 

2 The supercritical regime 10 

2.1 Modified equilibrium problem 10 

2.2 Determination of the Sj's 15 

2.3 Steepest descent analysis (supercritical case) 17 

2.3.1 The outer parametrix 18 

2.3.2 The local parametrix near a* 19 

2.4 Error analysis in the supercritical case 21 

2.5 The supercritical kernel and proof of Theorem 1.1 26 



1 



3 The subcritical regime 27 

3.1 Modified equilibrium problem (subcritical case) 28 

3.2 Steepest descent analysis (subcritical case) 29 

3.2.1 The local parametrix near b* 30 

3.3 The subcritical error analysis 31 

3.4 The subcritical kernel and proof of Theorem 1.2 35 

A The detailed analysis of Theorem 2.1 36 

A.l Background material 36 

A. 2 Complete proof of Theorem 2.1 for the case K = 37 



1 Introduction 

The physical motivation behind studying Hcrmitian random matrix ensembles is as a model of the Hamilto- 
nian for complex systems, where the eigenvalues of the random Hcrmitian matrix represent the energy levels 
of a system without time-reversal invariance [31]. 

Let A be a fixed Hcrmitian matrix. Wc consider the set of all n x n Hcrmitian matrices M endowed with 
the probability measure 

^(dM) = J_ e -« Tr ^( M )- AM )dM; Z n := [ e~ n ^ v ^- AM UM, (1-1) 

Z n J 

where rfM is the entry-wise Lebesgue measure and the integration is over all Hcrmitian matrices. 

When A = (no external source) and V(M) = M 2 /2, (1-1) describes the Gaussian Unitary Ensemble, 
or GUE. For A ^ and V(M) = M 2 /2, this measure arises in the study of Hamiltonians that can be 
written as the sum of a random matrix and a deterministic source matrix [16]. We are specifically interested 
in small-rank sources of the form 

A = diag(a, . . . , a, 0, . . . , 0) (1-2) 

r n—r 

assuming that either r = 0(n 7 ), < 7 < 1 or r is finite (in which case we define 7 := 0). The ratio of r to 
n, which is asymptotically small, will be denoted as 

«:=-. (1-3) 
n 

Peche [29] studied the limiting distribution of the largest eigenvalue in the Gaussian case (V(M) = M 2 /2) 
under these assumptions and found three distinct behaviors. In the supercritical case, r eigenvalues are 
expected to exit the bulk and are found to distribute as the eigenvalues of an r x r GUE matrix. For the 
subcritical case, the largest eigenvalue is expected to lie at the right band endpoint and behave as the largest 
eigenvalue of an n x n GUE matrix. In the critical case, when the outliers lie at the band endpoint, the 
distribution for the largest eigenvalue is an extension of the standard GUE Tracy- Widom function [30] which 
arises when r = (see also [1, 4, 6]). 
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One of the primary goals of random matrix theory is to determine universality classes of matrix ensembles, 
that is, find different probability measures on the space of matrices for which the spectral properties are 
the same in the large- n limit. With this goal in mind we consider more general functions V(M), with 
our specific assumptions listed in Section 1.3. Basically, we assume V(M) is a generic single-gap analytic 
potential with sufficient growth at infinity. We show that Peche's results [29] hold for these general potentials 
in the supercritical and subcritical cases. The universality of the critical case will be considered elsewhere 
[8]. 

1.1 Definition of the supercritical, subcritical, and critical regimes 

Let g(z) be the ^-function associated with the orthogonal polynomials with potential V(z) (see, for instance, 
[21] or [23]). It may be written as 

g(z) := / log(z - s)p min (s)ds, (1-4) 



where p m m is the unique probability measure minimizing the functional 

?\p] : = / V{s)p{s)ds- { { p(s)p(s')log\s- s'\dsds'. (1-5) 



We will assume this equilibrium measure p m i n is supported on a single band \a,0\ (see Assumption 1.1 (v) 
in Section 1.3). Define 

Pi(z) := -V(z) + 2g(z) + h, (1-6) 

P 2 (z) := -V{z) + az + g(z) + h, (1-7) 

P 3 (z) := -Pi(z) + P 2 (z) = az - g(z) - h + h- (1-8) 

Here g{z) and l\ are uniquely determined by the conditions that P\(z)± is purely imaginary on the support 
of the equilibrium measure and has negative real part on its complement in K (here the subscripts ± denote 
the boundary values from above/below the real axis). How l 2 is chosen will be described at the end of this 
section. 

It is also known that 3ig(z) is a continuous function on M and harmonic on the complement of the 
support of p m in (up to a sign it is also known as the logarithmic potential in potential theory). 

Definition 1.1. Define a c to be the (unique) value of a so that P 2 {P) = 0. 

The uniqueness is promptly seen because P 2 {0) = —V'(j3) + a + g'(/3); in fact the effective potential Pi 
is known [22] to satisfy 

P{(z) = 0(z - & (1-9) 
and in particular P{{(3) = and hence P 2 ((3) = a — g'(f3) = a — ^V'(f3). Thus the critical value of a is given 

by 

a c =g'(P) = -V'W). (1-10) 

We first have 
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Lemma 1.1. The critical value a c = g 1 (ft) is positive. Moreover g'(a) < 0. 

Proof. From (1-4) we see that g'(z) = J p mm (s)ds is positive for z > ft. It is also known that the 
density p m ; n vanishes like a square root at the endpoints a, ft and hence the integral representation of g 1 (ft) 
is convergent and immediately shows it to be positive. Similarly g'(a) is negative. Note that this proof 
does not require the support to consists of a single band as long as we understand ft = sup supp p m i n and 
a = inf supp/9 min . □ 

Lemma 1.1 implies that there is no loss of generality in studying only the case a > since there is always 
a positive critical a c (and a negative one); the negative case (a < 0) is equivalent to the positive case by 
replacing a H> —a, V(z) H> V(—z). 

Lemma 1.2. The critical point structure o/SRP^z) is: 

• For a > a c , -ft P$(z) is strictly increasing on M. \ [a, ft]; 

• For a = a c , ^RP^(z) is strictly increasing onR\ [a, ft] and ^RP^(ft) = 0; 

• For < a < a c , ^tP^(z) has unique local minimum onR\ (a, ft). This minimum occurs at a point 

b*e(ft,oo). 

Proof. From the representation (1-4) of g one sees immediately that 

g"(z) = - f r -t—5P mi n(s)ds (1-11) 

Jm. \ z ~ s ) 

which shows clearly that for z G K \ supp p m i n the real part of g is concave downward. Thus ^Rg'(z) is 
decreasing in M \ supp p mm ; moreover, from 

9'{z) = f -^—p mi n(s)ds (1-12) 
Jr (z - s) 

we see that -Kg' is negative for z < inf supp p m ; n = a and positive for z > sup supp p m i n = ft- 
From the definition we see that P^(z) = a — g'(z) and hence we infer: 

• For a > a c — g'(ft) > 0, ^.(P^(z)) = a - g'(z) > g'(ft) - g'(z) is positive on [ft, oo) therefore 5ftP 3 
is strictly increasing. On the other hand a — g'{z) is clearly positive on (— oo,a] because a > and 
-g' > from (1-12); 

• For a = a c , 5ftPg(/3) = and 5ftPg(z) is a monotonically increasing positive function on (ft,oo), and 
a monotonically increasing positive function on (— oo,a]. Therefore there is a single critical point of 
?ftPs(z) at z = ft. As 5ft Pg'(z) = —5ft g" (z) > this must be a minimum; 

• For < a < a c , ^RP^(ft) < and P^(z) — >• a > for z — > oo; moreover 5RPg(z) is a monotonically 
increasing function on [ft, oo), and a (monotonically increasing) positive function on (-co, a]. Since 
P^{ft) < there must be a unique point 6* > ft where P^b*) = 0. As 5RP^'(z) = -$tg"(z) > this 
must be the local minimum (or, equivalently, the global minimum on (ft, oo)). 
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□ 

Wc can now define four regimes: supercritical, subcritical, critical, and jumping outliers. We define the 
sub-critical and critical regimes first. 

Definition 1.2. The matrix model specified by (1-1) is in the subcritical regime if a < a c and P2(%) < 
P 3 (b*) for all x>b*. 

Definition 1.3. The matrix model specified by (1-1) is in the critical regime if a = a c and P 2 (x) < P2(P) 
for all x > (3. 

Now the supercritical regime can be efficiently defined as the remaining cases, except the small — 
codimension one — cases that wc distinguish as the "jumping outlier regime." 

Definition 1.4. The model is in the supercritical regime if Pi has a unique point of global maximum on 
{x > max{/3, b*}} at a point x = a* G ffi and any of the three conditions below is satisfied: 

• a> a c . 

• a = a c and Pi{P) < P2{x) for some x > (3. 

• < a < a c and P 3 (b*) < Pii^) for some x > 6*. 

Note that a* is always greater than (3 and b* . 

If the global maximum of P2 on (max{/3, 6*}, 00) is attained at several distinct points then we will say 
that we are in the jumping outlier regime that also includes the following remaining case. 

• < a < a c and Pzix) = Ps(b*) for some x > b* . (The case x = b* cannot occur for regular V .) 

In the present paper we consider the supercritical and subcritical regimes. The critical and jumping 
outlier regimes will be considered elsewhere [8]. 

The definition of the supercritical regime is complicated and the reader may wonder whether the above 
definitions ever hold in actual examples. It is however not difficult to engineer a situation where they do 
occur, explained in the following example 

Example 1.1 (Second and third bullet in Definition 1.4). Consider a potential V such that a new spectral 
band (i.e. interval of support of p m in) *s about to emerge. Then P\ has a local maximum —E at xq > /3 
outside of the main band which is slightly negative but small in absolute value. It is simple to arrange 
examples where E is arbitrarily small. Since a — a c we have 5i P 3 (/3) = (see (1-10) J and since is 
concave upwards, it must be increasing for x > (3. On the other hand P2 = P\ + P3 and thus 

(P 2 (x ) - P 2 (/3)) =-E + ®(P 3 (x ) - P 3 (/?)). (1-13) 

Since E can be chosen arbitrarily small, and since ^(i-^xo) — Pa(/3)) > 0, we see that necessarily we can 
have the situation described in the second bullet. By a continuity argument on a < a c , this also provides an 
example for the third bullet since P 3 (b*) < P 3 {!3) ~ P2W) < p2{%o)- 
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Though we do not consider in this paper, one can also create an example in the jumping outlier regime. 

We show in Proposition 1.1 that if V is convex, then we are either in the super or subcritical depending 
on a > a c or < a < a c , respectively; in particular, in this case, the situation described in the second and 
third bullet points of Definition 1.4 cannot occur. 

Proposition 1.1. Suppose that V(z) satisfies Assumptions 1.1 and in addition it is convex (V" > 0). Then 
(i) for a > a c the model is supercritical and the maximum of Pi at a* is nondegenerate; 
(ii) for < a < a c the model is subcritical; 
(iii) there is no jumping outlier regime. 

Proof. It is known that the convexity of V is a sufficient condition for the support of the equilibrium measure 
to be a single band [a,0\. Moreover from Lemma 1.1 we see V{[3) = 2g'(/3) > > V'(a) = 2g'(a). Note 
that Pj are all real-valued in [/3,co) (the cut of the logarithm runs in (— oo,/3]). 

(i) We then observe that 

P'i' = -V" + g" < (1-14) 

since both — V (by assumption) and g (by (1-11)) are concave downward. Pi may have at most a single 
global (nondegenerate) maximum in [/3, oo) because P'^ifi) = a — ^V'((3) = a — a c > 0. This also proves (iii). 

(ii) If < a < a c , Pi strictly decreases on [/J, oo) because P^iP) = a — i-V'(f3) = a — a c < 0.; Also we have 
P 2 (b*) = P 1 (b*) + P 3 (b*) < P 3 (b*). Therefore P 2 {x) < P 3 (b*) for all x > b* as in Definition 1.2. □ 

It should be noted here that, contrary to the work done for V = z 2 /2, the position of a relative to a c 
is not sufficient (for general V) to define the critical and subcritical regimes. If, however, V is convex (for 
example an even monomial with positive coefficient) then by Proposition 1.1 the position of a relative to a c 
determines the supercritical/subcritical regime as in [14, 3, 15, 1, 29]. The secondary conditions in Definition 
1.4 are dealing with whether the Lagrange multiplier li in the effective potentials Pj can be chosen such 
that the off diagonal entries of the jump matrices for the deformed Ricmann-Hilbcrt problems that we will 
construct in Sections 2, 3, and 4, decay to zero at an exponential rate. The problem of finding necessary and 
sufficient conditions on V(x) and a for the matrix model to be in the supercritical/subcritical regime is quite 
difficult, as much as it is difficult to find necessary and sufficient conditions for V(x) to be a single-band 
potential. We now specify the constant li. 

Definition 1.5. The constant li in Definition 1-7 will be chosen as follows: 

• In the supercritical case, the constant li is chosen so that the unique global maximum of Pi{z) on 
(/3,oo) is zero (i.e. $tP 2 {a*) =0). 

• In the subcritical case, the constant 

h:=-h+h (1-15) 

is chosen so that P 3 (b*) = 0. 
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1.2 The kernel and its connection to multiple orthogonal polynomials 



Let p m (Ai, . . . , A m ) be the probability density that the n X n matrix M chosen using (1-1) has eigenvalues 
{Ai, . . . , A,„} (here m < n). Then, the m-point correlation function is i? m (Ai, . . . , A m ) := ( n ^' m y p m (Xi, ■ ■ ■ , A m ). 
Brezin and Hikami [16, 17, 18, 19] showed that, in the Gaussian case, the m-point correlation functions can 
all be expressed in terms of a single kernel K(x, y): 



i? m (Ai,...,A m ) =det(K(Xi,Xj))i 



(1-16) 



Zinn-Justin [32, 33] extended this result to the case of more general V(M). Bleher and Kuijlaars [13] rewrote 
the kernel in terms of multiple orthogonal polynomials, a significant result because it allows one to analyze 
the asymptotic behavior of these polynomials via the associated Riemann-Hilbert problem. 

This approach was followed by Aptekarev, Bleher, and Kuijlaars [14, 3, 15] in the Gaussian case when 
the matrix A has two eigenvalues ±a, each of multiplicity n/2. When a is sufficiently large the eigenvalues 
of M accumulate on two disjoint intervals (the supercritical case). As a decreases the two bands collide (the 
critical case). Below this critical value of a, the eigenvalues accumulate on a single interval (the subcritical 
case). Related behavior also appears in the theory of nonintersecting one-dimensional Brownian motions; 
see, for instance, Adler, Orantin, and van Moerbckc [2] for the critical case. 

In general, the existence and number of bands on which eigenvalues accumulate for large-rank sources for 
general V(M.) is a complicated problem. For more on this question see McLaughlin [28] in which the quartic 
case V(M) = M 4 /4 is worked out. Bleher, Delvaux, and Kuijlaars [12] have studied the external source 
problem with two eigenvalues of equal multiplicity and where V(M.) is a sum of even-degree monomials 
with positive coefficients. The external source with a finite number of different eigenvalues with various 
multiplicity for supercritical case has been considered in [25]. 

The starting point of our analysis is the Riemann-Hilbert problem associated to the multiple orthogonal 
polynomials. Suppose Y(z) is a 3 x 3 matrix-valued function of the complex variable z satisfying 

'Y(z) is analytic for z ^ R, 
(\ 




Y + (x)=Y_(x) 



Y(z) 




n(V{x)—axy 



o 
l 






for x € 



(1-17) 



Here Y±(x) := lim^o Y(a;± ie) denote the non-tangential limits of Y(z) as z approaches the real axis from 
the upper and lower half-planes. Whenever posing a Riemann-Hilbert problem we assume (unless otherwise 
stated) that the solution has continuous boundary values along the jump contour when approached from 
either side. Under our assumption (iv) in Section 1.3, the unique solution Y(z) can be written explicitly in 
terms of multiple orthogonal polynomials of the second kind (see [14], Section 2). In the case of two distinct 
eigenvalues a and 0, the kernel K n (x,y) may be written in terms of the function Y(z) as 



K n {x,y) 



-l n (V{x)+V{y)) 

2iri(x — y) 



(0 1 e na y) Y{y)- x Y{x) | 

,0, 



(1-18) 
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In the technical analysis of this Riemann-Hilbert problem we use and improve certain ideas introduced by 
Bertola and Lee [9] to study the first finitely many eigenvalues in the birth of a new spectral band for the 
random Hermitian matrix model without source. 

We note here that Baik [7] has recently expressed the kernel K n (x,y) in terms of the standard (not 
multiple) orthogonal polynomials. This offers an alternative method for approaching the problem we consider 
here. Based on this approach, Baik and Dong [5] have obtained the universality result similar to ours, for 
the case of finite r but possibly for non-degenerate eigenvalues, i.e. A = diag(ai, a2, ■ ■ • , a r , 0, • • • ,0). Since 
the rank of the matrices involved in this alternative formulation grow with r, analyzing the Riemann-Hilbert 
problem (1-17) seems more feasible if r is allowed to grow sublincarly with n. 

1.3 Assumptions on A and V(z) and results 

First we gather the assumptions we will make in the rest of the paper. 
Assumption 1.1. We make the following requirements 
(i) a > 0. 

(ii) A is a small-rank external source of the form (1-2) with either r a fixed positive integer or r = 0(n 7 ), 
with < 7 < 1. When r is fixed we say 7 = 0. 

(Hi) V(z) is real-analytic. 

I iv Inn - — ; — = 00 and lim - — ; — = 00. 

|z|-voo l0g(l + Z 2 ) \z\^oo log(l + z 2 ) 

(v) V(z) is a single-band potential (for example it can be convex). 

(vi) The density of the equilibrium measure ofV(z) has square root decay at its two endpoints (i.e. it is 
regular in the sense of [23]). 

(vii) For the supercritical regime, P%(z) behaves quadratically near a* . Specifically, 

-V'(z) + a + g'{z) = -c(z - a*) + 0{{z - a*) 2 ) as z a* (1-19) 
for some constant c > 0. 

Assumption (i) is for convenience, as the case when a < is equivalent by sending a — > — a and V(z) — > 
V{—z). Regarding assumption (ii), in the general case when A has m > 2 distinct eigenvalues the kernel 
can be written in terms of multiple orthogonal polynomials associated to an (m + 1) x (rn + 1) Riemann- 
Hilbert problem, which is beyond the scope of this paper. Assumption (iii) allows us to use the nonlinear 
steepest-descent method for Riemann-Hilbert problems, while (iv) guarantees the existence of the multiple 
orthogonal polynomials needed to ensure the Riemann-Hilbert problem has a solution. 

Assumption (v) avoids the necessity of using Riemann-thcta functions for the solution of the outer model 
Riemann-Hilbert problem. We expect similar results to hold gencrically in the multi-band case. 
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Both (vi) and (vii) are genericity assumptions. Assumption (vi) allows us to use Airy parametrices near 
the band endpoints. Assumption (vii) produces Hermite (or Gaussian) behavior of the outlying zeros. Our 
results are computations of the large n behavior of the kernel function (1-18). We explicitly compute the 
kernel in a neighborhood of a* for the supercritical regime and in a neighborhood of b* for the subcritical 
regime. In the remaining portions of the complex plane, our result is that the kernel function converges to 
the kernel for the classical orthogonal polynomial problem with respect to V(x). In particular, our results 
include that, away from the a* and b* , the standard universality classes apply (i.e. a sine kernel in the bulk 
of the spectrum, and Airy kernels at the edges) . 



Theorem 1.1. Suppose V(z) and a satisfy assumptions (i)-(vii) and definition 1.4 of the supercritical 
regime. Let £ x and £. y be the local coordinates corresponding to x and y near a* as defined in (2-42). 
Uniformly for Cx,Cy * n compact sets we have the following asymptotics for r = Cn 1 for some C > 0, 
< 7 < 1. 

K n (x(C x ),y(Cy)) = e-^t^fftM^^-i/^ (k? ue & x , Q + ©(n^)/ 2 )) , (1-20) 
where Ps(x) is given by (1-8), f(z; n) is defined in (2-56), fc^-i * s defined by (2-62), k = r/n, and 

Kf UE (( x Xy) ~ ^(C* - Co)g^i(Cy - Co) - H { r r \(( x - Co)^ r) (Cy - Co) ^(C-CoJ'-SfC-Co) 2 (1 _ 21) 

Cx — Cv 

is the kernel for r eigenvalues of the Gaussian Unitary Ensemble of scale r centered at £o> which is defined 
by the change of variables (2-58). Here H^ r \C) are the reseated monic Hermite polynomials satisfying the 
orthogonality condition (2-62). 



The presence of exp(— nP^{x)/2) in (1-20) does not affect spectral properties of the kernel (because 
it amounts to a conjugation of the kernel by a diagonal operator) and therefore the implication is that 
asymptotically the eigenvalues near a* are equivalent to those of a scaled r x r GUE problem; if V(x) is a 
quadratic potential this agrees with the results of [29]. 

Remark 1.1. If the critical point of P2 at a* is more degenerate, Pv{z) = 0((z — a*) 2k ), then one may 
follow similar steps as in [9] and [11] and conclude that the relevant statistics of the outliers are determined 
by the kernel of a unitary ensemble with potential given by a polynomial of degree 2k instead of a Gaussian, 
namely, obtained from a deformation of the Freud orthogonal polynomials. 

Theorem 1.1 shows that, as expected, the equilibrium measure has no mass near a* . We have the mean 
density of states 

Pn (x(( x )) = lim -K n (x(C x ),y(Cy)) - ^P« 1/2 U r) (C x ) + 0{n-^' 2 r- 1 )) (1-22) 

'V— 1 
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where pr (£) is the mean density of eigenvalues for the r x r GUE ensemble. If r is fixed, this quantity is 
C(n -1 / 2 ) for large n, and if r = n 1 , using that pf\C,x) - > ^- 1/4 — £| as r -> oo, we find 

lim p„(x(C,)) = ^^^V^Cf. (1-23) 

In either case, our conclusion is that for large n the mean density of eigenvalues is asymptotically small (of 
order k 1 / 2 ) in the neighborhood of a* chosen in the theorem. See Chapters 5 and 6 of [21], Theorem 1.1 of 
[14], and Theorem 8.1 of [28] for similar results regarding the derivation of the asymptotic mean distribution 
of eigenvalues from the kernel. 



Theorem 1.2. Suppose the pair (V(z),a) satisfies Definition 1.2 of the subcritical regime. There exists a 
closed disk of fixed radius centered at b* such that, for x and y in this disk, for large n, and for r = Cri 1 for 
some C > and < 7 < 1, there is a c > such that 

K n (x, y) = 0(n-^/ 2 e- cn ). (1-24) 



Acknowledgments. The authors would like to thank Jinho Baik, Ken McLaughlin, Sandrine Peche and 
Dong Wang for several illuminating discussions. Wc thank Baik and Wang for sharing their unpublished re- 
sults. M. Bertola was supported by NSERC. R. Buckingham was supported by the Taft Research Foundation. 
V. Pierce was supported by NSF grant DMS-0806219. 

2 The supercritical regime 
2.1 Modified equilibrium problem 

In this section we will use a positive integer K; the general statements are valid for any K but we will choose 
(for future use) 

K > max jy— !-,oj (2-1) 

where, we recall, 7 is the exponent of growth of r = Cn 1 for some C and < 7 < 1. 

We will need to build a perturbation of the equilibrium problem that leads to the definition of g(z); we 
will denote by q the resulting g-function of this perturbation scheme. The construction, rather involved, 
will be broken down in steps. 

The unperturbed equilibrium measure is supported on the single interval [a, (3] (by assumption) with 
external field V(z). Recall that a* is lying outside of [a, 0\ and we fix a compact interval J containing [a, (3] 
in its interior and a* J. 
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Proposition 2.1. For any K £ N there is a neighborhood of the origin in (n,S) £ C 1+K such that the 
equilibrium measure a(x)dx of total mass 1 — k for the external field 

K 5 

V(z) := V{z) - 6V(z), 6V(z) := k ln(z - a*) + « £ _ J (2-2) 

zs supported on a single interval [a(n, S), f3(n, 5)] still contained in the interior of J: the endpoints a(n, S), /3(k, 5) 
are analytic functions of the specified variables. 

Furthermore the (normalized) g-function of this problem 

q{z) := / ln(z-w)P^-dw (2-3) 

J I - K 

converges uniformly over closed subsets not containing [a, 0\ to the unperturbed g-function. 

Remark 2.1. In this proposition we treat the deformation parameters k,S as independent from each other; 
later on, in Proposition 2.2, they will be uniquely determined in terms of the sole parameter k. 

Proof. It is well known (see, for example, [22]) that 

• q(z) is analytic for z ^ (— oo,/3], where /? = sup supp(cr); 

• q(z) has continuous boundary values and satisfies 

g+(z) - Q-(z) = 27ri, z <E (-co, 67), 

6 ' ~ (2-4) 

(l-K){ 0+ (z;K)+Q-{z;K)) = V(z)-Klog(z-a*)-y2 T , — - h, ze {a,[3) 

— ' 2 (z — a*y 

j=i ' 

where a = inf supp(5) with the real axis oriented left to right; 

• g(z; k) = log(z) + O ^-^ as z — s- oo. 

Vice versa, the g-function may be characterized by the scalar Riemann-Hilbert problem (2-4) with the 
additional requirement that 5g+ is a nondecreasing function. 

To show the analytic dependence of q on k, 5 we proceed as follows: define the function 

R(z):=((z-a)(z-f3))^ 2 , (2-5) 

where the principal branch of the square root is chosen so R(z; n) = z + 0(1) as z — > oo. Taking the 
derivative of (2-4) with respect to z and using the Plcmclj formula gives 

R(z) ^V / ( S )-n/( S -a*) + j:f =1 jS 3 /(2(s-a*)^) 
( z ) = T~71 \ / 1 \ u i \ ds 2 " 6 

2lTl(l — K) J 5 (S — Z)K+(S) 

where R+(s) refers to the limit in s from the upper half-plane. The large- z expansion of (2-6) along with 
the condition 

/ 1 \ 

(2-7) 



li 



^ives two conditions on a(/c), /3(k): 

V'(s) K/js - a*) + - a*Y +1 ) 

R+(s) 



ds = (2- 



and 

1 /•^'( S )- K /(.s-a*) + E-= 1 J^/(2( S -a^+ 1 ) 

sets = 1 — /t. (2-9) 



2m J 5 R.+ (s) 

These two equations uniquely determine 3, /3 as analytic functions of the parameters by the implicit function 
theorem. 

The inequality > remains valid, using a continuity argument, for suitably small values of k, 5 

because it is valid (with the strict inequality) for the unperturbed ^-function (by our initial assumption). 

Therefore the expression (2-6) yields a bona fide g function for the modified external field V in a neigh- 
borhood of (k,S) = (0,0). 

The expression for q may be obtained by integration; specifically 



B(z) = f fl'(s)ds-4 

J a 



(2-10) 



and i\ is also determined by the requirement that g(z) = ln(z) + Q(z 1 ) (without the constant term). 
Explicitly 

ii = liin (J fl'(s)ds-lnz^ (2-11) 

which expression shows that l\ is also analytic in the parameters, given the already proven analyticity of g'. 

The statement about the convergence follows easily by noticing that q' converges to g'(z) as desired (note 
that they both have behavior 1/z + 0(z~ 2 )). This is best seen from the integral representations and the 
already proven analytic dependence on the deformation parameters. □ 

In parallel with the definition of the functions P, we shall define 



Vx{z;k) := -V(z) + 2(1-k)q(z;k)+£ 1 , 
T 2 (z;k) := -V(z) + az + SV + (1 - k)q(z;k) + l 2 , 
T 3 (z;k) := -T 1 {z;k) + V 2 {z;k) 



K 



V{z) := V(z)-SV(z), 6V{z) := k ln(z - a*) + k ^ 



2(z - a*y 



(2-12) 
(2-13) 
(2-14) 

(2-15) 



where l 2 as in Definition 1.5 and is independent of the deformation. 

Using deformation and continuity arguments we have that (for a sufficiently small value of the deformation 
parameters k, 8) the real part of V 2 has a global maximum in a neighborhood of a* . The main tool in the 
analysis of the supercritical case when r grows shall be the next theorem. 
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Theorem 2.1. There exists a conformal change of coordinate p = p{z; k, S) fixing z = a* (i.e. p(a*; k, 5) = 
) that depends analytically on the parameters k, 6 such that 

K s 

V 2 {z) := -V(z) + az+(l- k)q + h + 2k\u(z - a*) + V J —— (2-16) 

~{ ( z ~ a F 

can be written as 

K 



V 2 (z; k, 5) = -hp - a) 2 + 2 K \np + b + V ^ (2-17) 
2 L — ' o-> 



3 = 1 H 

where the parameters a = a(ft,<5), b = b(n, 5) and 7 = 7(k;<5) are analytic functions of the indicated 
parameters. Furthermore the Jacobian 

9 -Z (2-18) 
85 

is nonsingular in a neighborhood of the origin (i.e. for k and 5 sufficiently small). 

Proof. To simplify the notation we set a* = (up to a translation this entails no loss of generality). We can 
write Vi as 

k s 

V- 2 = -f(z; K, 5) + 2k ln(z) + V -4 . (2-19) 

3=1 

By the definition of a* (which is now translated to 0), the function /(z; k, S) has the property that 

f(z;0,0) = ^z 2 (l + O( Z )) ,C>0 (2-20) 

and hence 

f(z; k, 6) = ^-z 2 (l + 0{z)) + 0(k, 5). (2-21) 

Let us fix any (smooth) curve in the parameter space n(t), 5(t) and denote by a d its tangent vector; we then 
must show the identity (2-17) for t near 0. We suppress the notation of the dependence on k,S for brevity 
in what follows, with the understanding that /(z), a, jj, b all depend on them. 

In this part we only sketch the main idea, leaving a full proof for Appendix A. Let D(r) be the open disk 
of radius r > and let fli be the Banach manifold of univalent, analytic functions p : O(r) — > C which 
fix the origin p(0) = 0; this is a closed Banach submanifold of all univalent analytic functions because the 
evaluation map is continuous. Define now 

M :=fti x C K+1 = {p= (p,a,b,7), p G Oi, a, b G C} (2-22) 

which is naturally also an infinite dimensional Banach manifold. We are going to show that the ordinary 
differential equation in d that derives from (2-17) is intcgrablc on Ai; taking the implicit differentiation of 
(2-17) we obtain 



-df(z) + 20k ln(z) + Ef=i ^ = (a - P + 1% £f=i ^r) 8p + db- da{p a) + 2dn In p + £f =1 ^ 
=^> dp 



-df(z) + 2d K ln U) + Ef =1 ^ + db + da(p -a)- d -J 



a-P+f-Eli^ 



p ^j=i P '+ L 

-p K df(z) + 2 P K In U) 8k + Ef =1 ^ + P K (P ~ a)0a - Ef=i <>- ,!>" " 
9p = p ^ - . (2-23) 

ap K+l _ p K+2 + 2Kp K _ j ljP K -J 
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Formula (2-23) should be regarded as defining a vector field on M, and this flow together with p(z; 0, 0) = 
J 2f(z; 0, 0) gives an initial value problem. To see this we have to remember that the tangent space to M. 
consists of 

TM := O x C A '+ 2 (2-24) 

where Hq stands for the Banach vector space of bounded analytic functions on D(r) (without the requirement 
of being univalent) mapping to 0. The denominator vanishes generically at K + 2 values pf, since dp must 
be an analytic function, the numerator must vanish at the same points and this yields a linear system for 
the K + 2 values db,da, d~fx, . . . ,djK- To sec how this works in more detail, let pj be the roots of the 
denominator in (2-23) 



K 



K+2 



-K+2 



ap 



K+l 



31]P 



K-j _ 



For k, a, 7 sufficiently small all the roots pj belong to the disk 
there are corresponding values z\ , . . . , zk+2 ■ 

The linear system that determines da, c?7 is then 



-U(p-Pj)- (2-25) 

3=1 

(r) where p(z) is univalent and therefore 



Pf 



-df(zt) 



21n[^ 

Pi 



K K 



K 



3 = 1 Z l 



pf{ Pl - a)da + pfdb + J2 pf~ 3£> ^ = . 1=1,.-., K + 2. 



(2-26) 

What we want to see is that this system determines 9a, db, d'y as analytic functions of k, 0,7; to sec this we 
observe that the coefficient matrix of the linear system (2-26) is 



p¥(P2 - a) 
Pk +2 (pk+2 - a) 



Pl<+2 



Pf" 1 
Pf" 1 



PK+2 



H(z e ) :=pf 2 In 



da ' 
db 

d^K 

97x 
dn - df(z e ] 



H{zi) 
H{z 2 ) 



K+2 



K 



E"' 

3=1 



K 



dSi 



(2-27) 



(2-28) 



Solving this system by Cramer's rule yields da, d~j as symmetric functions of the roots pf, moreover it is seen 
that the determinant of the linear part is simply the Vandermonde determinant A.(p) := Ylj<e<K+ 2 (Pj — P^) 
and since the determinants in the numerators of Cramer's formula will also vanish whenever two roots 
coincide, it follows that the ratio is actually analytic on the diagonals pi = pk, I 7^ k. 

We have thus proved that da,d~f,db are analytic symmetric functions of the pgs; it is well known that 
the ring of analytic symmetric functions is generated by the elementary symmetric polynomials in the roots, 
namely, the coefficients of the polynomial (2-25). This means that da, dj, db are actually expressible in 
terms of analytic expressions of a, k, 7. 

In order to complete the proof we should check that the vector field determined by (2-23) is Lipshitz 
with respect to the Banach norm of TAi; the check is rather straightforward but lengthy and a detailed 
analysis is deferred to App. A in the simplified case K = 0. After this, the existence and uniqueness of the 
integrated flow follows from standard theorems in Banach spaces. 
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Jacobian at the origin. To compute the Jacobian at the origin (k, S) = (0, 0) we have to set 



V2/(z;0,0) , a = K = 5 J ; = 7j =b = . 
Taking now de to mean d$ e we find the equations 



0, 



(2-29) 



(2-30) 



Since we want p(0; K, S) = we must impose that dip in (2-30) vanishes at least of order z at z = 0; this 
yields a linear system for the coefficients di b, di a, de"/j and in particular 



^7,(0,0) - 0,j>£ 

9,7,(0,0) = 1 

d aj (0,0) = *,j<l 



85 



1 * -k 
1 * 



(2-31) 



with the * denoting some expression which is not relevant to us; the above Jacobian is then triangular with 
1 on the diagonal, and hence it is invcrtible in a neighborhood of k = 0, 6 = 0. □ 

2.2 Determination of the S^s 

We now introduce the rescaled variable (k > 0) 



C 



Co 



Let 



9h(C) ■= 4VC^4 + ln(C + Ve^) + £ + % 



2 2 
1 - 2 log 2 



(2-32) 

(2-33) 
(2-34) 



be the ^-function for the Gaussian Unitary Ensemble. It admits an asymptotic expansion of the form 



g H (Q :=lnC + ]T 



■HI 



(2-35) 



We define the constants Cj by the requirement 



K C (H) 



It is easily verified that the c, are polynomials in Ci 



(2-36) 



Proposition 2.2. 77ie parameters 8 are uniquely determined as Puiseux series of t/k by the requirement 



K 

W = - f (C - Co ) 2 + 2k ln(v^C) + « E 



- + b = --(C - Co) 2 + 2 K ln(V<) + b + 2k]T ^-(2-37) 



Moreover we have 



3 = 1 S 

Co = 0(\/«), b = <D( K ) , 



(2-38) 
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Proof. Recall that £o depends on both k, S analytically; although it is possible to give more detailed infor- 
mation about this dependence, it will not be necessary to the end of establishing the present proposition. 
We need to solve the nonlinear system 

K-i lj ( K ,$) = 2Kcf\( ) , j = l,...,K (2-39) 

for the unknowns 6. The local solvability of the system around 6 = is guaranteed if we can guarantee the 
nonsingularity of the Jacobian matrix. But this system can be rewritten 

7-2^)=0, D:=diag(| 2) 5,...^). (2-40) 
Since S- H ^ is analytic in S one promptly sees that the Jacobian is 

J:= d 4-2 K ^ (2-41) 
86 OS 

and hence det J = 1 + 0(k§). This guarantees that there is a polydisk |k| < C\, \\8\\ < C-i (for suitable 
constants) where the system admits a solution in Puiseux series (i.e. analytic in \/~k). It is also clear from 
(2-40) that 5 = 0{ni). Thus, Co( K 7<K K )) i s still of order 0(tJ~k) and b{K,8{n)) is still O(k), since all these 
depend analytically on 5, which is not of higher order than k. □ 

For future reference and definiteness we collect the result of the above discussion in the theorem below 



Theorem 2.2. There exists a conformal change of coordinate Q(z] k) of the form 

{{z; k) = = -±=C(z - o*)(l + 0(z - a*)) , C > (2-42) 

V K v K 

and a choice of 5 = for the deformed potential (2-2) in Puiseux series of ^JH such that 

k ^ c( H ~> 

^ 2 (r; K ,J( K )) = --(C-Co) 2 +2 K ln(^C) + b + 2 K ^— . (2-43) 

j'=i 

The functions £q(k), (3(k), 5(k) admit a Puiseux expansion and are of orders 

Co^O(V^), = O{k), 6 = 0(k). (2-44) 
The expressions cf 1 ^ are polynomials of degree j in £ determined by the formula (2-36) . 
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2.3 Steepest descent analysis (supercritical case) 

We make the following change of variables from Y(z) to W(z): 

Y(z) 















W(z;k) 


- 


i 









I o 








e -9V 






















e f (V-2az) 



(2-45) 



x < 


























& %(2V 2 -V l -2 V ) 



(2-46) 



= «;ln(z - a*) + kEJLi a(»-a«y 
:= k1iik + b + k£ h , (£# := -1 - 2 In 2) 



(2-47) 
(2-48) 



The constant (in z) r\ is chosen to carefully balance other constants later on in the study of the local 
paramctrix; we recall that b = O(k) has appeared in Theorem 2.2. which is also the source of the k\uk 
term. The constant £h was introduced in (2-34) and it is the Robin's constant for the equilibrium problem 
associated to the quadratic potential. See Figure 1 for a visual on the different regions f2j's. The exact 
choice of the outer lenses is given below in the proof of Lemma 2.1(d). The inner lenses are chosen in the 
standard way for the 2x2 Riemann-Hilbcrt problem for (non-multiple) orthogonal polynomials. 

The new matrix W(z) satisfies a new Riemann Hilbert Problem which can be directly evinced from the 
one for Y and is of the form W + (z) = W_(z)V( w )(2i) with jumps 



V (w) (z) 



'1 e n ~Pi( z ) e ™("P2(z)-»7) N 




z G 3fii n dfl 6 , 

z G (dQ 1 n dQ 2 ) U (9^5 n dQ 6 ), 
z G (<90 2 n 9Q 3 ) u (dn 4 n dQ 5 ), 

z G [a, (3], 

z g dn 2 n dfi 5 , 



(2-49) 
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Figure 1: The oriented jump contour T for W(z) and the regions fij in the supercritical case. 



and the asymptotic conditions 

( ~ ) , z -> oo 

/ e -nSV(z) o Q 

W(z) = (analytic) | 10 

o e nSV{z) I (2-51) 

as z — >• a* . 

The orientation of the contours is given in Figure 1. Here we have used the factorization 
'1 1 1\ /l 

1 = 1 — 1 J | 1 1 | ( -1 | ( 1 1 | [ 1 1 | . (2-52) 

,0 1/ \0 




2.3.1 The outer parametrix 



We will show below in Section 2.4 that the jump matrices for W(z) decay uniformly to constant jump 
matrices as n — > oo outside of small fixed neighborhoods of a, [3, and a* . These limiting jump matrices are 
the identity except on the band [a, 0\. We therefore define the outer parametrix *f?(z) to be the solution to 
the following Ricmann-Hilbcrt problem: 



*+(z) = *_(z) 



1 0^ 

-10 0) for z e (a,/3) 



= 1 + 



lj 

It is well known that the solution to this Ricmann-Hilbcrt problem is 



*(z) = U" 



-1/4 







— a 




1/4 



\ 




V 



u, u 



(2-53) 



(2-54) 
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where 

1/4 



U m ( '^JL ) = l (2-55) 



z-s-oo \ z — a 
and this function is cut along [a, /?]. 

2.3.2 The local parametrix near a* 

Special attention is needed near the point z = a* as the jump matrices do not decay uniformly to the identity 
near this point. According to the definition of a* as the point of maximum for P 2 and given that "P 2 is a 
deformation of P 2 we have 

:=/(*;«) 



P 2 (z; k) = - %{z - a*) 2 (l + 0(z - a*)) + 0(6) +2 K log(z - a*) + £ ^j. . (2-56) 
2 * — ' z — a P 

where C > and the deformation 0(<5) is some analytic function of z of the indicated order in n. Let D a * 
be a fixed-size circular disk centered at a* chosen small enough so that 



3? 



< |RPi(*)| (2-57) 



(recall that 9Wi > outside of the support of the equilibrium measure) inside the disk, and such that the 
disk docs not intersect the outer lenses. Orient <9B a * clockwise. 

We now apply Theorem 2.2: let the k constants , j = 1, . . . , k be specified by (2-36): in the local 
scaling coordinate £ we have 

r - 
nV 2 = --(C-Co) 2 + 2rlnC + 2r^] — +rhxK + nb. (2-58) 

j'=i ^ 

where the scaling coordinate £ has the following behavior on the boundary of the disk D a * 

C = C(n (1 ~ 7)/2 ) when z e <9!V. (2-59) 

We also recall (Theorem 2.2) that b = O(k) and hence nb = 0(r). This suggests the following definition for 
the model Ricmann Hilbcrt Problem of the local parametrix. 

Definition 2.1. The local parametrix within the disk D a * shall be the unique solution R(z) to the following 
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model Riemann-Hilbert: 



(l C 2r expf-§(C-Co) 2 + ^ + 2rEf =1 4?")\ 



R+(C) = R-(C) 



R_(C) 



1 

\0 

/ 1 e 71 ^ 2 '^ 

1 

\0 1 



( i ) 



R(£) = (analytic) 



C r exp(rEf=iV 



(2-60) 



as ( — S- 0. 
?7 := k1iik+ b-«4» (2-61) 



The behavior at C = (z = a*) is dictated by (2-51). We point out that the problem is essentially 2x2; 
moreover it will be shown below that it is a slight modification of the Fokas-Its-Kitaev Riemann-Hilbert 
problem for Hcrmitc orthogonal polynomials (sec [27] and [24], Section 3). 

Proposition 2.3. Let the rescaled Hermite polynomials h!> (£) be the family of monic polynomials satisfying 
the orthogonality condition 



where the ' are normalization constants. Then the solution to (2-60) is 



(2-62) 



K (H) 



R(C) = exp (--^Ais) H 13 (C)r rAl3 exp -l H — A 



k-13 



(2-63) 



where 



Hi 3 (C) 







2?ri 



OO 



^ r) ( S -Co)e-5 s2 \ 
—as 



2?ri 



^(C-Co) o 



V-l 



oo S-C 



oo trv. 



1 f°° H^s-^e 



-ds 



and 



Aia 



10 


l0 -1; 



(2-64) 



(2-65) 



The proof is a direct manipulation and it is left to the reader. 
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Now the well known asymptotics of Hcrmite polynomials can be written as a joint asymptotic expansion 
for £ and r large as follows 



Hi 3 (C) = e 



_ ^Ih A13 



I + O 



1 



kici + i; 



u 



C-Co-2 
C-Co+2 









\ 



U e rgH(C-Co)A ia -^e H Aia ^2-66) 



4/ C-C0+2 

C-Co-2 / 



1 i 
1 

-1 i 



"2-67) 



and gniO 1S given by (2-33). For large Q, we have the expansion of gn as in (2-36). The error term of 
0(l/r(|C| + 1)) in (2-66) is from the standard Airy paramctrix used in the Ricmann-Hilbcrt problem for 
Hermite orthogonal polynomials at Q = ±2. For more details on this calculation see equation (7.72) in [21] 
or equation (4.16) and Appendix B in [24], noting that the variable ( is rescaled by a constant factor. Then 
for large ( (such as on dB a * ) we can estimate 



1^ 



I 



oil 



exp 



K C (H) 

to(c-Co)-iogc-E^r 



1-13 



(2-68) 



O 



»K+1 



2.4 Error analysis in the supercritical case 

Let TD a and Bp be small, fixed, closed disks centered at a and (3 that are bounded away from the outer lenses. 



Orient the boundaries dB a and dBp clockwise. For z 6 D Q , let P^(z) be the Airy parametrix satisfying 



• P 



(z) has the same jumps as W(z) for z € 



. Pg (*)*(*)- 



O ( - ) for z e 

n 



The construction of the Airy parametrix is standard, involving Airy functions and a local change of variables. 
See [14] Section 7, for example, for an Airy parametrix for a 3 x 3 Riemann-Hilbcrt problem. The Airy 
parametrix P^( z ) is defined analogously for zeDp. 
We now define the global paramctrix ^°°(z) by 



The error matrix E(z) is given by 



*(z)R(C(z)), zeD,,, 

z eB a 



«■»<«>■ 



(2-69) 



z e Bp. 



E(z) := W(z)*°°(z) 



(2-70) 



Let r denote the contours given by the boundaries of the regions Qj in Figure 1. The error matrix satisfies 
a Ricmann-Hilbert problem with the following jumps: 



21 




Figure 2: The jump contours r( E ) for the Riemann-Hilbcrt problem for E(z) in the supercritical case. 



• For z outside the disks D a , V)p, and D a *, and excluding the band [a,/3]: 

V ( E )( Z ) = *(z)V( w )(z)*(z)- 1 , z e r n (B Q U Dp u d q *) c n [a,f3] c 

where V^ w \z) is given by (2-49). 

• For z on the boundaries of the disks <9B Q , 9D^, and <9B a *: 



'*(z)R(0*(«)-S 
V( E >(z)=<|pg(z)*(z)-\ 



zedB a *, 
z e <9B Q , 



• For z inside the disk D a * : 




z e r n D a * 



Furthermore, V^ E ' (z) = I on the contours 

[a,/3]n(© Q u© (3 ) c , rn: 



and T n 



(2-71) 



(2-72) 



(2-73) 



The jump contours r( E ) for E(z) are shown in Figure 2. 

We now show that all of the jump matrices in (2-71)-(2-73) are uniformly close to the identity as n — > oo. 

For the error bounds it will be convenient to split r( E ) into a compact component and a noncompact 
component 



-(E) < 



: u 5% u do a . u (rniv), 



p(E) ._ r (E)\ r 
1 N ■— 1 \ A 



E) 

C 



(2-74) 



We now gather the results we will need on the functions Pi(z), Pi(z), and P${z) defined by (1-6)— (1- 
Lemma 2.1. In the supercritical regime, the inner and outer lenses can be chosen so that: 
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(a) On the inner lenses outside of the disks around a and f3: The real part of P±(z) is positive and bounded 
away from zero for z G [(<9Q 2 H dil 3 ) U (dfl 4 n <9Q 5 )] n (B„ U By3) c . 

(7>^ On the real axis outside of \a,0\ and the disks around a and j3: The real part of P\{z) is negative and 
bounded away from zero for z G [(<9Qi H dtt e ) U (8Q 2 H 90s)] n (B a U %) c . 

fcj On £/ie rea/ axis outside of the band [a, /?] and a fixed distance away from a, ft, and a* : The real part of 
P 2 (z) is negative and bounded away from zero for z G [(<9^2 H Sfis) n (B Q U B^ U B a *) c ] U (<9f2i n df^)- 

('dj On t/ie outer lenses: For k sufficiently small, the real part of Pz{z) is negative and bounded away from 
zero for z G (dfii n <9ft 2 ) U (<9r2 5 n <90 6 )- 

Proof. Statements (a) and (b) follow from the analysis of the Riemann-Hilbert problem for the standard 
(not multiple) orthogonal polynomials (see, for instance, [21]). Statement (c) follows from the definitions of 
the supercritical region and a* . 

For (d), we begin by choosing the outer lenses used to define W(z). Fix n = 0. Note that 

P 3 (ft) = -i\08) + P 2 (/3) = P 2 (J3) < P 2 {a*) = 0. (2-75) 

The second equality uses the fact that 3? Pi is zero on [a,0\ and Pi is real for x > f3\ the inequality follows 
since a* is the location of the global maximum of P 2 {z); and the final equality is true by the choice of the 
Lagrange multiplier l 2 . Thus, there is a fixed radius neighborhood of j3 in which RePz(z) < for real z. We 
choose the outer lenses to be a circle centered below a whose right-most endpoint passes through the real 
axis at some point on (f3, a*). We choose the circle big enough such that %IP 2 is negative on the real axis to 
the left of the circle. This is always possible due to Assumption 1.1 (iv). 

We now show that the outer lenses are descent lines of ReP3(z). Clearly the real part of az decreases 
as we move to the left along the lenses. Note that Keg(z) = / log|2 — s\p m i n (s)ds where p m \ n {s) is the 
associated equilibrium measure. Now for any s G (a, j3), log \z — s\ is increasing as z moves to the left along 
the lens (one can see this clearly by drawing a circle that is centered at s and is tangent to the outer lenses 
at the right-most point). Since p m i n (s) is positive for s G (a, ft), Kcg(z) increases as z moves to the left 
along the lenses. This shows (d) for n = 0. Since S(/t) = a + O(k) and /3 = /3 + O(n), (d) also holds for k 
sufficiently small. □ 

We now present the results we will need for the modified functions Vi(z; k), V 2 (z; k), and Vz{z; k) defined 
by (2-12)-(2-14). 

Lemma 2.2. For K sufficiently small (refer to Figure 1 and Figure 2): 

(a) On the inner lenses outside of the disks around a and f3: The real part of V\{z; k) is positive and 
bounded away from zero for z G [(dfl 2 n dfl 3 ) U (dfl 4 n dO, 5 )] n (B Q U Hp) c . 

(b) On the real axis outside of[a,/3] and a fixed distance away from a, j3, and a* : The real part ofV\{z;n) 
is negative and bounded away from zero for z G fl dQ$) U (dQ 2 D dQ$)] D (B Q U B^ U B a *) c . 
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(c) On the real axis outside of the outer lenses and a fixed distance away from a* : The real part ofPi{z\ k) 
is negative and bounded away from zero for z € (<90i n dfle) l~l (D Q *) C . 

(d) On the outer lenses: The real part ofV 3 (z; k) is negative and bounded away from zero for z g (dtti f~l 
dQ 2 ) U (dQ 5 n dQ e )- 

(e) For real z inside D a * ; Let gaiO be defined by (2-33). Then the real part of 

K (H) 

Pi (*;«) + Kg H (0 - k log C - k ^ ^j- 

is negative and bounded away from zero for z € T D D a * . 

Proof. Parts (a) through (d) follow from Lemma 2.1(a)-(d) together with the convergence g g guaranteed 
in Proposition 2.1, the boundedness of log(z — a*) and (z — a*)'^ , j = 1,...,K outside of D a * , and Sj (0) = 0. 
For part (e), first note that comparing the two expressions (2-13) and (2-58) for V-2,{z;k) gives 

K \og(z - a*) - «logC + ^ 2 (z -a*V ~ K ^^T = 2 /(z; k) _ 2 & ~ Kln ^~ 4K " ^) 2 - (2-76) 
j=i 1 ' 3=1 ^ 

By the choice of D a « (see (2-57)), for k sufficiently small we have 

"l , b 



2/( z ; K ) - - KlnVK- -(C- Co) 2 



< \RPi(z)\ for z S D a *- (2-77) 



Write 



Pi (z; k) + Kg H (C) - K log C - k ~kf 

i=i ^ 

A 5- i,c (ff) (2-78) 

= -V + 2(1 - k)q + h + ng H + n\og{z - a*) ~ n\og( + V J -— - k V ^— . 

=Pi+o( K ) > ^ J 



Has real part bounded above by |5?Pi| 

For k small, by the convergence q — > 5 in Proposition 2.1, the first group of terms on the left-hand side, 
— V + 2(1 — k)q + £1 + Kgn, is within O(k) of Pi, which has strictly negative real part for 2 G D a *. Since 
the real part of the second group of terms is strictly less than the magnitude of the real part of Pi , part (e) 
follows. □ 

We are now in a position to bound the jumps 'V^(z) of the error problem. 

Lemma 2.3. In the supercritical regime, for large n, 

(a) Outside the disks D a , Hip, and D a * ; There is a constant c > such that 

V( E )(z;«)=I + 0(e" c "), zeffl 
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(b ) On the boundary of D a * : 

V( E )(z; k) = I + O (n-^-^/ 2 ) + O (n^-^ K+ V) , z e dB a * . 

(cj On the boundaries o/D Q andD^; 

V (e) (z;k) = 1 + (n" 1 ) , ze9D Q U5D^. 

fdj Inside D a * ■' 

v (E) (z) = i + 0( e - c ™), zernB a *. 

Proof. Part (a) follows from (2-71), Lemma 2.2(a)-(d), and the boundedness of Part (b) follows from 

(2-59), (2-68), and the boundedness of &(z). Part (c) comes from the construction of the parametrices 
P^(z) and Pf}{z) (see, for instance, [24]). 

For part (d) we consider the jumps (2-73) inside the disk D a *. Looking at the formula (2-63) for R(C), 
it appears there may be a problem at £ = 0. However, note that 

'1 (*)l2 0> 

V (e) (z;k) =*(z)e-5^ Al3 Hi3_(C)e- r(to(c) ~ £ f-) Al3 I 1 ojx 



,0 1 (2-79) 



, n p r(g H {0 

wherein 



^5 Al3 Hi3-(0 _1 e 5<HAls *(«) _1 , ^ ernD„., 



/ K c {H) \ 

{*)i2 = expUVi{ziK)+rg H (0-rlogC-rJ2 J ^-j ■ ( 2 " 80 ) 

Now (2-79) together with Lemma 2.2(e) and the boundedness of ^(z) inside D a * establishes (d). □ 

We can now show that the error matrix E(z) is uniformly close to the identity. 
Lemma 2.4. In the supercritical regime, for n large, 

E(z) = l + o(n~^-^/ 2 \ 

uniformly in z. 

Proof. From Lemma 2.3(b)-(d), 

VW(z) = l + 0(n-^l 2 )+0(ni-^ K+ V), zeT^. (2-81) 

The first error term always dominates or matches the second term if the nonnegative integer K is chosen so 
(2-1) is satisfied. Then, for n sufficiently large there exists a constant c such that 

HV (E) -I|l i2(r C E )) + HV( E > -HU^)) < cn-^. (2-82) 

Also, from Lemma 2.3(a), for n sufficiently large there is a constant c such that 

l|V(E) ^(rCf ) + H V(E) - * C6 ~ Cn > (2 " 83) 

The result follows by a standard technique that consists of writing the solution to the Riemann-Hilbert 
problem in terms of a Neumann series involving V^ E ' — I (sec, for instance, [24] Section 7.2 or [26] Section 
3.5). □ 
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2.5 The supercritical kernel and proof of Theorem 1.1 

Proof of Theorem 1.1. Recall that the kernel is defined by (1-18): 

-L 2 n{V(x) + V{y)) 

Kn{x ' v) = 2*i( X -y) ([nyr'Y^^+e-y [Y(y)-'Y(x)} 31 ) . (2-84) 

We consider local coordinates ( x and Cy m ©a* ■ While the function Y(z) has a jump in this region, the first 
column of Y(z) does not (see the Riemann-Hilbert problem (1-17)). Therefore we can pick x and y to be in 
a convenient region. We choose x and y to be in f2i as defined in Figure 1. 
From the transformation (2-45), we see that 

[Yiyy'Yix)] 21 = [W^)" 1 W(z)] 21 exp (n ((1 - K )g(y) + (1 - n)g(x) + SV(x) +£{)), (2-85) 

[Y(y)- 1 Y(x)] 31 = [W(y)- 1 W(x)] sl exp(n((l-K)Q(x) + 6V(x) + 6V(y) + l 2 -v)) (2-86) 

for x and y in Oi . We have 

W(z) = (1 + 0(n- (1 ^ )/2 )*°°(z) = (1 + 0( ? i- (1 - 7)/2 )*(z)R(C(z)) (2-87) 

for z G D a * . From (2-42), we have 

^(j/)- 1 *^) = 1 + O ((&, - C,)« 1/2 ) • (2-88) 

We define the functions Qi(z; re) to be: 

Qi(z;k) := -V(z)+2(l-/s)fl(«;/s) + ^i(/c), (2-89) 
Q 2 (z;k) := -T/(z)+az + (l-K)g(z;K) + ?2, (2-90) 
Q 3 (z;k) := az-(l-«)fl(jB;«)-*i(/s) + f 2 . (2-91) 

Now combining (2-85), (2-87), (2-88), (2-63), detR(C) = 1, and noting the 0{k) error terms from (2-88) arc 
subsumed by the ©(n - ^ 1-7 ^ 2 ) error terms from (2-87) gives 

e -f(v(x)+v(v)) [Y(y)-i Y (x)] 21 = (O ((C - C*> 1/2 ) ■ ^ r) (C ~ Co) + 0(n-( 1 ^)/ 2 )) e »W, (2-92) 

where 

1 1 K c {H) 

(*) + ^Qi(y) + - KiogC, - 

1 1 3=1 ^ (2-93) 

=\Px{x) + \Pi{y) + 0{K\riK). 

The last equality is shown by noticing that rearranging the terms in (2-43) we have 

k (J?) 

<JV(x) - K log C - k "V = o ~ ax -t 1 - K ^ ~ ^2 - ^ (C - Co) 2 + ^ In K + b) (2-94) 

3=1 C * 

Since the real part of Pi(z) is negative for z near a* (Lemma 2.1(b)), for k sufficiently small the real part 
of the exponent in (2-92) is negative. 
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Define 

P r GUE (C„ C„) := ^ (H^\Q - ( Q )Hl; r \(t y Co) - H { ;\{C, X ~ - Co)) ■ (2-95) 
From (2-86), 

e -^{v{ x )+v{y))+nay [ Y {y^Y {x)\ ^ = (F r GUE (£ x , ( y ) + O (n~ (1 ~ 7)/2 )) e n <**>, (2-96) 
where from (2-86) 

V( ) V( ) K c {H) K c {H) 

2 2 ^ CI V f=i CI (2-97) 
+ (1 - k)q(x) + k£ h . 

Using now (2-94) (for both x and y) and rearranging the terms we find 

1 r k 
(**)= o (! " K )fl(^) ~ o(Cy - Co) 2 - ^2 + + b 



(1 - k)q{x) -ax--(( x - Co) 2 - h + tlnK + bj + k£ h - i] (2-98) 



= --(& - Co) 2 - -(£, - Co) 2 - ifisCx) + ifi3(y). (2-99) 

where we have used the definition of 77 := ft In ft + b + kHh (2-48). From this and the fact that (2-92) is 
exponentially decaying in n shows 

Kn{x ' y) = er ^Z + - Q y) V) (^ GUE ^'^) + °( n " (1_7)/2 )) er^-^ 2 -^-^ 2 . (2-100) 

One final application of (2-42) (to switch x — y to £ x — Cy) an d Proposition 2.1 (to show convergence of Q3 
to P 3 ) gives (1-20). 

□ 



3 The subcritical regime 

We now take V(x) and a so Definition 1.2 of the subcritical regime is satisfied; that is a < a c and the function 
5ft P2 has no global maximum on R\ [a, /?]. In this case 5ft P3 has a (unique) global minimum at z = b* (with 
value zero, as per our choice of I2 in Definition 1.5). We will show that almost surely there are no outliers. 
We will freely reuse the same notation from Section 2.1 for new objects which played a similar role in the 
analysis of the supercritical regime. To begin, fix 7 € [0,1) and again set K to be the smallest nonnegative 
integer satisfying 

K > maxjy— !-,oj (3-1) 
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3.1 Modified equilibrium problem (subcritical case) 

The procedure here parallels closely the one followed in the supercritical case, and hence we will only state 
the results since their proof does not differ significantly from the other case. 

Let J be a closed subset not containing the point b* and containing [a, j3] in its interior; we recall that 
b* (a) > (3 for < a < a c . 

Proposition 3.1. For any K £ N there is a neighborhood of the origin in (k,S) £ C 1+k such that the 
equilibrium measure a(x)dx of unit total mass for the external field 

K S 

V(z):=V(z) + SV(z), SV(z) := nln(z - b*) + ac]T 2(z ' (3-2) 

3=1 

is supported on a single interval [oc(k, 6), /?(«, 5)} still contained in the interior of J: the endpoints a{n, 5), P(k, 5) 
are analytic functions of the specified variables. Furthermore the g -function of this problem 

q(z) := J ln(z — w)a(w)dw (3-3) 

converges uniformly over closed subsets not containing [a, (3] to the unperturbed g-function. 

The proof is identical to that of Proposition 2.1: the only difference is that now the modified equilibrium 
measure is of unit total mass, rather than of mass 1 — k. We next re-define the three functions "Pj's; 
the definition is subtly different from the previous (2-12), (2-13), (2-14) and hence there is a possibility of 
confusion for the reader. The advantage is that we will be able to recycle many of the previous computations. 











M*) 


:= -V(z)+2g(z)+h, 




(3-4) 


V 2 {z) 


:= V 1 (z) + P 3 (z) = -V(z) + 6V(z) + az 


+ Q(z) + l 2 


(3-5) 


Vz{z) 


:= az-Q{z) + 25V(z)+h-£ 1 =V 2 (z)- 


-Vx{z) 


(3-6) 










V(z) 


:= V(z) + 6V{z), 5V{z):=K\og{z- 


>*>+e 2(2 !v,^ 

3=1 


(3-7) 



It is important to point out the change of sign in the definition of V, relative to the supercritical case. 
We also remind that (i,Q are analytic functions of k, 5, while l 2 is the constant mandated in Definition 1.5. 

Theorem 3.1. There exists a conformal change of coordinate p = p{z;n,8) fixing z = b* (p(b*]K,S) = 0) 
that depends analytically on the parameters K, S such that 

K . 

V 3 (z) :=az-Q(z; K ,5) + l 2 -£ 1 + 2 K \n(z-b*) + y i J —- (3-8) 

*-rf (z — b*y 
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Figure 3: The regions f2j and the oriented contour T for the subcritical case. 



can be written as 

1 K 
V 3 ( z;K J) = -(p-a) 2 + 2K\np + b + Y,^ (3-9) 

where the parameters a = a(K,o), b = b(n,d) and 7 = t*(k; 6) are analytic functions of the indicated 
parameters. Furthermore the Jacobian 



(3-10) 



07 
85 

is nonsingular in a neighborhood of the origin (for k sufficiently small). 
Theorem 3.2. There exists a conformal change of coordinate Q(z; k) of the form 

C(z; k) = ^fi = ~^C{z - &*)(1 + 0{z - b*)) ,C>0 (3-11) 
and a choice of 8 = 5(k) for the deformed potential (2-2) in Puiseux series of \fn such that 

K C (H) 

V 3 (z; k, S(k)) = -- (C - Co) 2 + 2 K ln(^C) + b + ^- . (3-12) 

The functions Co( K )j S(k) admit a Puiseux expansion and are of orders 

Co = 0(Vk), P = 0(k), 6 = 0{k). (3-13) 
The expressions c~- are polynomials of degree j in Co determined by the formula (2-36). 

3.2 Steepest descent analysis (subcritical case) 

The regions fii, i = 1, ...6 are defined in Figure 3. Following the example of the supercritical case, we 
introduce the same new matrix W as in (2-45). Note, however that the definition of the regions fi/s now 
follows Figure 3, the constant I2 follows Definition 1.5 in the subcritical case and SV is given in (3-2) instead. 
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The new matrix W(z; k) satisfies the same jump conditions (2-49) as in the supercritical case (with, of 
course, the new definitions of the Vj(z\ k)'s and fij's) as well as the same asymptotic condition (2-50). Due 
to the new definitions of Vj's (3-4), (3-5), (3-6) , the behavior near z = b* is now different: 



'10 

W = (analytic) ( c - nSV W 

i0 e ^v{ z ) 

as z -> b*. 



(3-14) 



The outer parametrix problem is the same as in the supercritical case, so we again define the outer parametrix 
solution \l>(z;re) as in (2-53). 



3.2.1 The local parametrix near b* 

Define ©6* to be a fixed-size circular disk centered at b* which is small enough so that it does not intersect 
the inner lenses, and 

5R(P 2 )<0 for zeD k .. (3-15) 

This last condition is possible for k sufficiently small because -due to the definition of li in Definition 1.5 for 
the subcritical case- P 2 (&*) < 0. We now use Theorem 3.2; in a local coordinate centered at b*, the analytic 
part of Vs(z; k) behaves the same way (quadratically with a maximum at the origin) along the imaginary 
axis as Vi (z; k) did along the real axis in the supercritical regime. The scaling of C is analogous to (2-59); 



C = O (n* 1 " 



•7)/2 



when z £ dl 



(3-16) 



k (H) 



P 3 = f(z; k) + 25V (z) + l 2 - l x = -|(C - Co) 2 + 2^ log C + 2k ^ ^- + K In k + b 



(3-17) 



Definition 3.1. The local parametrix within the disk D&* shall be the unique solution R(z) to the following 
model Riemann-Hilbert problem: 



(\ 



R+(C) = R-(C) 



= R-(C) 



o 



( \ ) 



1 -C 2r cxp^-§(C-Co) 2 +^ + 2rE; =1 ^r 

\o o i 

/l \ 

1 -e nV3 , (el, 

V° 1 / 



(3-18) 



R(C) = (analytic) 



/I 

r r exp(-r£) =1 ^ 



7 exp ( r X^Li ^jhr j j 



as C ->• 0. 
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This problem is almost the same as (2-60). Analogously to (2-63), the solution is (£h 
(2-34)) 



k JH) 



R(0 = ex P (-^Ajs) h 23 (o c rA23 ^ll^-'-E^r] A2 3 



1 — 2 In 2 as in 



(3-19) 



where 



H 23 (C) := 



f 1 




V 



o 

H«(C-Co) 



2?ri 



^i(C-Co) 



-1 
1 





^ r) (s-Co)e"^ 



s-C 



8-C 



-ds 



-ds 



and A 
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/ 



(3-20) 

Again the polynomials Hm\() and the normalization constants fcm are defined by (2-62). The analysis in 
the supercritical regime leading to (2-68) applies here as well, leading to 



( \ 



O 



2fe+2 



(3-21) 



3.3 The subcritical error analysis 



Let B Q and B^ be small, closed disks of fixed radii centered at a and (3 that are bounded away from the outer 
lenses and B/,* . Orient the boundaries &D a and dHp clockwise. Let P Vj and P^9 be the Airy parametrices 
constructed in B Q and B^, respectively (see Section 2.4). Define the global parametrix ^f°°(z) by 



*°°(z) := < 



*(2)R(C(«)), zeBf, 
P<g(z), zeB Q , 



The error matrix E(z) is given by 



E(z) := W(z)* 00 (z)" 1 . 



(3-22) 



(3-23) 



Let r denote the contours given by the boundaries of the regions Qj in Figure 3. The error matrix satisfies 
a Riemann-Hilbert problem with jump matrix V( E ) (z) on the contours r( E ) shown in Figure 4. The form 
of the jump matrix is as follows: 

• For z outside the disks D„, H)p, and Of,*, and excluding the band [a, /3]: 



v (E) (z) = *(z)v (w) (z)*(z)" 1 , ze rn(D Q uD^uD 6 .) c n [a,p] c 



(3-24) 



where V< w )(z) is given by the formulas in (2-49). 

• For z on the boundaries of the disks <9B Q , 9B^, and <9B;,*: 



'*(z)R(C)*(z)- 1 , zedB b ., 
V (E \z) = \ P ( A a . ) (z)*(z)" 1 , z G dB a , 



(3-25) 



z G 51 
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Figure 4: The jump contours r( E ) for the Ricmanii-Hilbcrt problem for E(z) in the subcritical case. 



For z inside the disk 



V (E \z) 



*(z)R(C) 



/l 





i"Pi(z) 



1 I RiC)- 1 *^)- 1 

\0 1 

1 e nVl{z) 0\ 
|0 1 R(C)- 1 *^)- 1 , 
1/ 



z e (3fiin90 6 )nD(,, 
z g {dft 2 n 9fi 5 ) n D b * 



(3-26) 



• Furthermore, V^ E \z) = I on the contours 

[a, 0\ n (D Q u D/3) c , rnD„, rno^, (afiin9fi 2 )nD t ., and (<9ft 5 ndfi 6 )niv. 

We now show that all of the jump matrices in (3-24)-(3-26) are uniformly close to the identity as n — >• oo. 
Here are the results we will need for Pi, P2, and P3, analogous to Lemma 2.1: 

Lemma 3.1. 7n t/ie subcritical regime, the inner and outer lenses can be chosen such that 

(a) On the inner lenses outside of the disks around a and f3: The real part of P\{z) is positive and bounded 
away from zero for z G [{dtt 2 H dn 3 ) U (<9ft 4 H <9Q 5 )] n (B Q U ©,g) c . 

^ On i/ie reaZ axis outside of [a, /3] and the disks around a and j3: The real part of P\(z) is negative and 
bounded away from zero for z E [(<9Q X n 9fi 6 ) U {dVL 2 n 90s)] n (B a U B'/}) c . 

(cj On the outer lenses outside of the disk around b* : For k sufficiently small, the real part of P 3 {z) is 
negative and bounded away from zero for z G [(9f2i PI d£l 2 ) U (dfl^ H <9f2g)] H . 

(c() On i/ie real axis outside of the outer lenses or on the real axis inside Bf,*: For k sufficiently small, the 
real part of P 2 {z) is negative and bounded away from zero for z G (9Sli l~l d£l§) U [(dfl 2 H 3^5) H O5*]. 

Proof. Parts (a) and (b) follow from the analysis of the Ricmann-Hilbcrt problem for the standard (non- 
multiple) orthogonal polynomials (for example, [21]). 
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For (c), first note Ps(b*) = 0. Following the proof of Lemma 2.1(d), the outer lenses (defined in this 
regime to be a circle centered below a and passing through b*, that is big enough such that Sfti-2 is negative 
on the real axis to the left of the circle) are descent lines of Rcpj(z) for k sufficiently small. The result 
follows. 

For (d), start with k = 0. Consider real z to the left of the outer lenses. From Lemma 3.1(c), RePj < 
at the left- most point of the outer lenses. Thus KePs(z) < for such z since Rep3 is a strictly increasing 
function for z G (—00, a). Since RePi(z) is also negative here by construction, this means ReP2(z) = 
Rc(Pi(z) + P^{z)) is also negative. This is also true for real z inside HV by (3-15). Next, consider real z to 
the right of the outer lenses. By Definition 1.2 we have Rep2(z) < RePs(b*) = 0. Along with the fact that 
Re Pi (z) < here shows the desired result . □ 

Next come the necessary results for Vi(z; k), Vi{z\ k), and V^z] k) defined by (3-4)-(3-6). This lemma 
is analogous to Lemma 2.2 for the supercritical regime. 

Lemma 3.2. For k sufficiently small: 

(a) On the two inner lenses outside of the disks around a and {3: The real part of V\(z\ k) is positive and 
bounded away from zero for z G [(e?fi 2 H dfl 3 ) U (9Q 4 n dQ 5 )} n (B a U B fj ) c . 

(b) On the real axis outside of [a, /3] and the disks around a, (3, and b* : The real part ofV\{z; k) is negative 
and bounded away from zero for z G [(Sfii n dfl e ) U (dfl 2 n dfl 5 )} n (D Q UD^U B b *) c . 

(c) On the outer lenses outside o/Df : The real part ofVs^z; n) is negative and bounded away from zero 

for z g [(0fti n dn 2 ) u (dn 5 n dn 6 )] n Dg*. 

(d) On the real axis outside of the outer lenses and D5* ; The real part ofViiz^ k) is negative and bounded 
away from zero for z G (dQ,\ n dfle) PI B^;* . 

(e) On the real axis inside TS>b* : The real parts of 



are negative and bounded away from zero for z G [(0f2i fl 8Qq) U (8^2 H 9^s)] H D^* . 

Proof. Parts (a) and (b) come from Proposition 3.1 along with the fact that log(z — b*) and (z — b*)~i , 
j = l,..., 2k, is bounded outside of O;,* . Part (c) comes from combining Lemma 3.1(c) with Proposition 3.1 
and the boundedness of log (2; — b*) and (z — b*)~i , j = 1, . . . , 2k. 

Part (d) follows from Lemma 3.1(d), Proposition 3.1 , and the boundedness of log(z — b*) and (z — fe*) -3 , 




k AH) 




k 1H) 



and 




□ 



These results allow us to now bound the jumps V< E > of the error problem. Divide r' E ' into a compact 
component T c and a noncompact component T N : 
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Lemma 3.3. In the subcritical regime, for large n: 

(a) Outside the disks D Q , D^, and : There is a constant c > such that 

V { V\z-K)=l + 0{e- cn ) 1 zeT™. 

(b ) On the boundary of HV : 

V< E ) (z; k) = I + O (n-^' 2 ) + O (V fe - 1+ ( fe+2 ^) , z e dO b > 

(c) On the boundaries o/D Q andUip: 

V (e >(z;k) = 1 + fij , zedB a UdB p . 

(d) Inside D&* : There is a constant c > such that 



v ( - E \z-,K) = i + o(e- cn ), zern: 



Proof. Part (a) is the result of Lemma 3.2(a)-(d) and the boundcdncss of ^f(z). Part (b) is from (3-16), 
(3-21), and the boundedness of *S?(z). Part (c) is from the construction of the parametrices P^' (z) and 
Pju(z) (see, for instance, [24]). 

For part (d), consider the jumps (3-26). By (3-19) for R(C), 

, A (*)i2 (*)is\ 

V^(z;k) =*(z)e-5 fffA23 H 23 _(C)e- r(!?H(c) -^ L)A23 1 O x 

\0 1 J (3-28) 

x e '"(ffH(0-^)A2 3H23 _^)-i e |^A 2 3^^- ) -i ) z G (gOjn^jnDfc,, 



and 



/I (*)ia 0\ 

V^(z;k) =*(0)e-5 fffA23 H 23 -(Oe- r(9ff(c) "^ L)A23 1 O x 

\0 1/ (3-29) 

x e r ( to («-^) A23 H 23 _(C)" 1 e5^ A23 *(z)- 1 , z e (9fi 2 nffi 5 )nD fc ., 



wherein 



Co- 



(*)i2 =exp [ riPi(z;n) - rg H (() +rlogC + r^-^j- 



(*)i3 = exp ( nV 2 (z; k) +rg H {() - r log C - r -^j- ] • 



(3-30) 



This along with Lemma 3.2(e) and the boundedness of *&(z) in IDV gives the result (d). □ 

We can now show that E(z) is asymptotically close to the identity. The proof of the following lemma 
follows that of Lemma 2.4: 

Lemma 3.4. In the subcritical regime, for n large, 

E(z) =I + o(n- (1 "T'/ 2 ) 

uniformly in z. 
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3.4 The subcritical kernel and proof of Theorem 1.2 

Proof of Theorem 1.2. Once again, recall the kernel (1-18): 

-±n{V(x) + V(y)) 

Kn{x ' y) = 2*i(x - y) ([ Y <*) + enaV [ Y ^) Y W]3i) ■ ( 3 " 31 ) 

While the function Y(z) has a jump for z e IV , the first column of Y(z) does not (see the Riemann-Hilbcrt 
problem (1-17)); observing the Riemann Hilbert Problem for Y -1 we also note that the second and third 
rows of Y -1 are entire functions. Therefore we can pick x and y to be in a convenient region. We choose x 
and y to be in f2i as defined in Figure 3. 

From the transformation (2-45) (which is the same in the subcritical case as noted at the beginning of 
Section 3.2) and using the new definitions (3-4), (3-5), (3-6), we see that 

[Y(y) _1 Y(x)] 21 = [W(y)- x W(x)] 21 cxp (n ( B (y) + g(x)-SV(y) + h)) , (3-32) 



[Yfo)- 1 Y(a:)] 31 = [Wiy^Wix)] gl exp (n ( S (x) + SV(y) + l 2 - r,)) 
for x and y in fii. As in the supercritical case, we have 

and 

*(y)- 1 *(x)=I + 0((C x -( y )K 1/2 ). 
Define Qi(z; n) to be Vi(z; k) without the logarithm or pole terms: 



Qi(z;k) 



-V(z) + 2q(z;k)+£i, 
—V(z) +az + q(z; k) + l 2 , 
az — q(z; k) + l 2 - l\. 



Recall (3-20) that 



R(0 = O(l)e V / asC^O. 



Now combining (3-32), (3-35), (3-39) gives 
and 

e -f (V(»)+v( v ))+no» [Y(y)- 1 Y(a;)] 31 = C((x - 2 /)n^( 1 '^/ 2 )e™ ( ** ) . 
Rearranging the terms from (3-8) and (3-12) we find 

K (H) 

5V(z) - k\u( -Kj2^7- = ^(9(z)-az -?-(<;- (o) 2 + Kln K +b-l 2 +l 1 ) 

.7=1 



(3-33) 

(3-34) 
(3-35) 

(3-36) 
(3-37) 
(3-38) 

(3-39) 

(3-40) 
(3-41) 

(3-42) 
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Using (3-42) we can rewrite 

(*) = \QiW + 562(1/) + ~(C V ~ Co) 2 ~ «lnv^- \ = \Px{x) + l -P 2 {y) + 0(«ln«) (3-43) 

and 

(**) = ifii(x) + is 2 (y) - ^(C y - Co) 2 - \ + kIh = \Pi(x) + \p 2 {y) + 0(«ln«). (3-44) 

Here we have used Proposition 3.1 to convert g to <?. Since 5ft-Pi(6*) < and also #tP2(b*) < 0, the theorem 
follows. □ 

A The detailed analysis of Theorem 2.1 

In this appendix we prove the existence of the local change of variables used in the supercritical and subcritical 
cases. We also demonstrate how the change of variables can be computed explicitly termwise. 

A.l Background material 

We start recalling that the space %{p{r)) of holomorphic functions on an open connected domain U> (a disk 
of radius r for simplicity) is a Banach space with respect to the sup norm. 

The theorem of existence and uniqueness for ODEs can be extended to any Banach space Ai . A sufficient 
condition for the intcgrability is the Lipshitz property, namely that we are given a (time-dependent) vector 
field V : A4 x J — > TA4 which is jointly continuous and Lipshitz. Let 

fti := {C : B(r) -»■ C, C(0) = 0, ||C||oo < oo, C univalent} (A-l) 

and 

fi := {C : 0(r) ^ C, C(0) - 0, ||C||oo < oo}. (A-2) 

Lemma A.l. The evaluation map of the inverse p^ 1 at a point is locally Lipshitz on More precisely: 
VCo € fii BC, P ,S>0 s.t. VC € C, |f | < p VCi, C2 G B s (Co) C fii (A-3) 

ICi _1 (0 - < ClICi - C2II00 = C sup |Ci(z) - Ca(*)| • (A-4) 

zGD(r) 

Proof. Note that fii is an open subset of the Banach vector space of bounded analytic functions on H>(r) 
containing the identity map. Therefore the Banach ball of radius S > centered at £o G Oi lies all within 
Sli for sufficiently small 5. 

First we note that the forward map is locally Lipshitz; that is, let zq € D(r/2), then any of the functionals 
(( n \z ) are Lipshitz 

itf°(*)-d B) (*>)i = 



2i?r 



(&(*)- C2(*))d* 



z|=2/3r (-2 — ^o) 



n+1 



< 



2?r 



,n+l 



IC1-C2 



(A-5) 
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Let now Co(-z) € f2i be univalent and let < 3p := infui =r /2 ICo(-z)|- Let £ be such that |£| < p. 
Let C11C2 be two maps in a ball around £0 of radius p (WQ — Coll < P) and consider (all integrations are 
on \z\ = \r) 



\z,-z \:=\Q\Z)-Q 1 (t)\ = 



< 



2tt 



2i7rjT Ci(*)-£ ziKJbW-Z 
*[(C««)(6-0-G(*)(Ci-0]d* 



(Ci - 0(C 2 - 



(A-6) 
(A-7) 

(Ci-o(c 2 -o ■ (A ~ 8 ' 

Since the derivative evaluation on the circle z = r/2 is uniformly Lipshitz, the above can be easily estimated 

by 



1 

2^ 



g IC(C 2 - CO + C( (C2 - Ci) + (C( - CjKi] d* 



1 

2^ 



g K(C 2 - CO + CjXCg - d) + (Cj - Ca)Ci] d* 

(Ci-0(Ca-0 



< 



cv 

inf |Ci-£| inf IC2 — ^1 

|«|=r/2 I I =r- /2 



IC1-C2 



(A-9) 



where C = sup{|£|, supi z |_ r / 2 |C(|j sw P\ z \=r/2 10 1}- Since Q are less than p away from we have that the 
two infima in the denominator are at least p, since |£| < p and 



, inf |0 1 > . inf „ I Co | - P = 3p - P = 2p. 

\z\—r/2 \z\—r/2 



(A-10) 
□ 



A. 2 Complete proof of Theorem 2.1 for the case K — 

Proof. In this case there are no i5's and no 7's; it should become clear that the general proof presents only 
notational complications, but is amenable to the same logic and hence the details are omitted. We also omit 
explicit reference to the dependence of f(z) on k for brevity. We want to have 

1 



-f(z) + 2/trnz = -|(p(z) - a) 2 + b + 2n\np(z) , 



(A-ll) 



where the goal is now to show that b = b(n), a = a(n) and p — p(z; k) are all analytic functions of k, with p 
being univalent in a neighborhood of z = and mapping the origin to the origin. Consider the differentiation 
of the above identity with respect to k: 



-f(z) +21nz = [ a- p+ -) p+ (p- a)a + /3 + 21np 



Solve for p and we find 



p(z;k) 



(p-q)d + b + 21n(f) +f(z) 



(A-12) 



(A-13) 



p — dp — K 

We want to view this equation as defining a vector field on a suitable Banach space that we define 
presently. Let fii be the Banach manifold of univalent, analytic functions p : D(r) — s- C which fix the 
origin p(0) = 0; this is a closed Banach submanifold of all univalent analytic functions because the evaluation 
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map is continuous. One only has to verify that if p n is a sequence of univalent analytic functions on D(r) 
converging in the sup-norm, the limit exists and it is still univalent. Define now 



M := fti x 



{p = (p,a, *>), (efli, a,be 



(A-14) 



Formula (A-13) defines a vector field on M: we will first explain it in coarse terms and then refine the 
details. 

The denominator to (A-13) has two roots pi(a, k), /^(fl, K ); since p must be also analytic, we must impose 
that the numerator vanishes at the same points, and hence 



a=A:= 



B 



det 


-M 


-2In(pi/*i) 1" 
-2ln(p 2 /z 2 ) 1 


det 


pi - a 1 
P2 - a 1 





det 


'(Pi 


-a) 
-a) 


-M)-2hi(p 1 /z 1 )" 
-/(z 2 )-21n(p 2 /z 2 )_ 


det 


pi - a 1 
p 2 - a 1 





Pi, 2 := 



a± Va 2 + 4k 



(A-15) 

(A-16) 
(A-17) 



Here 01,^2 are the counterimages of pi,p2, Zj '■= P~ 1 (Pj)- Note that the expressions have analytic 
continuations to the case p\ = p 2 : indeed they are symmetric functions of the roots and therefore they can 
be expressed in terms of analytic functions of a, n (which play the role of elementary symmetric polynomials 
in the roots). Therefore we consider the (time dependent) vector field V on the manifold M 



V([(p, a,b,K)])=( V (z),A,B) 



' (p-a)A + B + 21n(f)+/(z) ^ ^ 
p 2 - ap - k ' 



(A-18) 



It is to be pointed out that rj(z) — r/([p, a, b, k]; z) is a tangent vector to f2i, namely r]([p, a, b, k];0) = 0. 
The initial condition for the vector field is 



Po = (p(z; 0), a(0), b(0)) = v 7 ^), 0, 



(A-19) 

Therefore the proof shall follow if we show that the vector field V is integrable in some neighborhood of 
the initial point and for sufficiently small values of k and for this to hold it is sufficient to verify the Lipshitz 
property. 



Lipshitz property for V. To complete the proof it is sufficient to show that the vector field is locally 
Lipshitz in a Banach neighborhood of the initial condition. The initial p(z) is univalent in a small disk -say 
D(2r )- around z = because p'(0) ^ 0. 

By simple continuity arguments in the sup norm, there is a sup- neighborhood IA of p consisting of univalent 
functions on D(ro). 
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Wc therefore shall restrict a, k in such a way that \pj\ < rg; this guarantees that we can define the 
components of V. It is also quite clear that the restriction \pj(a, k)\ < ro contains a polydisk in a, k (here b 
is unrestricted). For example we can require 

H < To ■ W < Wo ■ ( A " 2 °) 

Thus, the neighborhood of the initial condition {po} x {0} that we will analyze is 

ZS~{(p,K) = (p,a,b,K)eMxR; P eU, \a\ < ^ , \k\ < ^| (A-21) 

The goal is now to prove that V is Lipshitz on 93 for any fixed t. The fact that A, B in (A-15, A-16) are 
Lipshitz functions follows from Lemma A.l. As for the first component, we have recall that the product of 
two Lipshitz bounded functions is Lipshitz, as well as the ratio if the denominator is bounded away from 
zero. By construction t](z) (A-18) is analytic and hence its sup-norm is achieved on the boundary of B(r) (by 
the maximum modulus theorem). By the restrictions we made on |a| and |«|, the denominator is bounded 
away from zero on 9D(r) and uniformly so with respect to the choice of p in the Banach neighborhood 14 of 
a. On the other hand, the numerator is clearly Lipshitz. □ 
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